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ABSTRACT 

Previous simulations of laser Doppler anemometer 
(LDA) systems have focused primarily on noise studies 
or biasing errors. Another possible source of error is 
oue to the choice of filter types and filter cutoff 
frequencies. 

In general, before it is applied to the counter 
portion of the signal processor, a Doppler burst is 
filtereo to remove the pedestal and to reduce noise in 
the frequency bands outside the region in which the 
signal occurs. Filtering, however, can Introduce errors 
into the measurement of the frequency of the input sig- 
nal which leads to Inaccurate results. 

With the results of this paper it is possible to 
evaluate errors caused by signal filtering in an LDA 
counter-processor data acquisition system and to choose 
filters for a specific application which will reduce 
these errors. 

NOMENCLATURE 

a s amplitude of Doppler burst peak, counts/sample 
period 

A,j c DC or background component of Doppler burst, 
counts/sample period 

F) lower cutoff frequency of the filter, cycles/ 
sample period 

F u upper cutoff frequency of the filter, cycles/ 
sample period 

fg Doppler frequency of the burst, cycles/sample 
period 

G maximum value of | H( f )| 

H(f) frequency response function of a given filter 
I H(f )| magnitude of H(f ) 

K comparison switch constant (generally set to 10) 


Kfl arbitrary constant (usually 3 or 4) for setting 
counter thresholds 

N number of cycles used in comparison testing 

N c number of cycles used in Doppler frequency 

measurement 

Ny number of fringes in the burst 

N bw noise bandwidth, cycles/sample period 

P n noise power 

P s mean signal power 

r comparison ratio, 1/2 or 5/8 

SNR peak signal to noise ratio 

T n time interval between zero-crossing time at the 
end of the first cycle to exceed TH2 and the 
zero-crossing time at the end of the last cycle 
in the comparison envelope 

Tfl C time interval between zero-crossing time at the 
end of the first cycle to exceed TH2 and the 
zero-crossing time at the end of the last cycle 
in the measurement envelope 

TK1 low counter threshold (K a x background noise 
level ) 

TH2 high counter threshold (5 x THl ) 

tg time of peak of Doppler burst 

V visibility of the burst 

x(t) time domain output of the photomultiplier 

<x(t)> expected value of x(t) 



o rms noise 

« phase angle of Doppler burst 

INTRODUCTION 

To date, most LDA data acquisition system analyses 
have focusseo on the errors caused by noise or biasing. 
This paper evaluates the errors caused by signal fil- 
tering through a computer simulation of the signal 
processing sequence. 

In general, a Doppler burst is filtered before 
being applied to the counter-processor to remove the DC 
component (pedestal) and to reduce noise from outside 
the frequency region containing the Doppler signal. The 
ideal filter for this application possesses three sig- 
nificant characteristics (1): linear phase (sometimes 

stated as constant time delay) in the passband, constant 
amplitude in the passband, and infinite attenuation 
outside the passband. While an ideal filter is theo- 
retically impossible, approximations to this ideal can 
be made. These approximations to the ideal are the 
general filter types analyzed herein. 

Filter-induced errors can be divided into three 
broad categories which correspond to the lack of per- 
fection in the previously mentioned characteristics. 
These error categories are as follows: skewing of the 

zero-crossings used in determining the input frequency, 
caused by nonlinear phase delay in the filter; biasing 
or preferential selection of certain frequencies in the 
passband, caused by ripples in the filter amplitude 
within the passband; and errors caused by insufficiently 
attenuatea noise from outside the passbana. 

Adrian and Earley ( 2 ) evaluated LDA counter- 
processor performance as a function of the noise attrib- 
utes of the signal received by the photomultiplier. 
Dopheide and Taux U) compared the counter-processor to 
a transient recorder but did note that "electronic fil- 
tering of the Doppler signal introduces the most impor- 
tant systematic error" to the data acquisition process. 
Hosel and Rodi ( 3 ) emphasized noise induced errors and 
determined that noise errors were a function of the 
signal to noise ratio (SNR), the noise spectrum band- 
width, and the counter control parameters. Except for 
Dopheide and Taux, when electronic filtering was men- 
tioned at all, it was assumed that the filtering being 
done was appropriate for the signals being measured. 

The emphasis of the work by Dopheide and Taux, however, 
centered on the comparison of two data acquisition 
techniques, rather than filter type and use. 

DIGITAL simulation of signal PROCESSING with the 
COUNTER-PROCESSOR 

Doppler Burst 

The output of the photomultiplier tube is assumed 
to be a Gaussian moouloteo sinusoidal wave of frequency 
equal to the Doppler frequency. It is further assumed 
that the seed particle velocity vector which generates 
this signal is parallel to the fringe normals of the 
measurement volume. The parameters the. describe the 
burst are the amplitude A s , the numbe • of fringes Nf 
(equal to the number of cycles between e~? intensity 
points), the visibility V, the Ooppler frequency fg, 
the time at the peak of the burst tg, and the phase 
♦ . The constant background level is A(j c . The 
photoelectron rate is then 


X(t) - A d( . ♦ A $ |l + V COS [2*fg(t - t # ) + ♦ ]} 

which is shown in Fig. 1(a) (in the time domain), and in 
Fig. 1(b) (in the frequency domain as represented by the 
magnitude along the positive frequency axis). Assuming 
a "short-term" counting condition (4), the number of 
counts (photon events) in a given interval ( t- ( dt/2 ) , 
t + (dt/2)) is approximately x(t) dt. The actual number 
of counts in each time interval dt has a Poisson 
probabi lity distribution. The probabi . ity of k counts 
in the interval dt is 

P(X,m) n e (2) 

where m is the expected count at time t 

m = x ( t ) dt (3) 

The Doppler burst of Fig. 1(a) with shot effect induced 
white noise added is shown in Figs. 1(c) and (d). 

Figures 1(e) and (f) she- the signal after filtering 
with a bandpass filter. Figure 1 illustrates the 
oesired elimination of the pedestal and the reduction 
of the noise. The number of visible fringes has a 
direct affect on tne apparent performance of the fil- 
ters. Figure 2 shows an exaggerated view of this 
affect. The Doppler frequency of the bursts was shifted 
upwards to the upper cutoff frequency of the filter to 
highlight the pedestal leakage. Figure 2(a) depicts the 
spectrum of a Dopple burst with Nf = 4, while Fig. 2(b) 
shows the spectrum of a Doppler Durst of identical fre- 
quency with Nf = 12. As can be seen in Figs. 2(c) 
and (d) (showing the signals of Figs. 2(a) and (b), 
after filtering), the pedestal in the 4 fringe case is 
much less attenuated than its 12 fringe counterpart. 

This pedestal leakage can result in frequency measure- 
ments which are much lower than the actual Doppler 
frequency. 

F i 1 ters 

The filters were simulated using classical analog 
filter oesigns, with the low pass prototypes being used 
to calculate the high pass filters using the usual 
transformation (5). The conversion to the digital 
uomain is done by sampling the complex frequency 
response of the analog filter which is equivalent to 
impulse invariance (6, p. 198). For this to be a valid 
proceoure the impulse response must decay to near zero 
in the time record. The bandoass filters were formed 
by cascading low and high pass filters of the required 
types. A naroware realization of these bandpass filters 
would require that an amplifier be inserted between the 
filters to maintain an acceptable signal level. 

Four five-pole filters were examined in an attempt 
to determine which of the three criteria were the most 
significant in terms of filter-induced errors. The 
Butterworth filter (Fig. 3(a)) was chosen because the 
low pass prototype provides a transfer function which 
is "maximally flat" in the passbana, although this 
characteristic is not maintained in the transformation 
to create the bandpass filter. The Bessel filter 
(Fig. 3(b)) low pass prototype exhibits a transfer func- 
tion which is "maximally flat" in terms of the time 
delay ’n the passband, but again, this highly linear 
phase characteristic is not maintained when the low pass 
filter is transformed to generate the bandpass filter 
(5). Finally, the Chebyshev filters were chosen as 



examples of filters with transfer functions of the low- 
est order which meet the specified attenuation criteria. 
The 1-dB filter (Fig. 3(c)) is an example of a filter 
with a fast fall off outside the passband, and the 0.1-dB 
filter (Fig. 3(d)) is an example of a filter with a low 
ripple. As a result of the low ripple, the Chebyshev 
0.1-dB filter also has a flat passcand magnitude. 

Signal Filtering Using Discrete Foirier Transform (OFT) 

For the digital simulation, the time record x(t) 
of the burst was represented as 1024 uniformly space 
samples. The Fourier transform of this record reveals 
the Doppler frequency. In this paper, the simulated 
frequencies are expressed in units of cycles per sample 
period. Consider a Doppler burst of frequency fp 
sampled 25 times per cycle with Nf = 8 visible 
fringes. The Doppler frequency of this signal is then 
1/25 or 0.04 cycles/sample period, and the burst width 
is Nr/' f o = 8/0.04 » 200 sample periods. 

A scaling factor is used to relate these numbers 
to actual signals obtainable from an '.DA sy .em. For 
example, with a scaling factor of 10® sample periods/ 
second, the 0.04 cycles/sample period signal would cor- 
respond to a signal of frequency 40 MHz. The burst 
width would then be 200 sample periods/10® sample 
periods/second or 200 nsec in a time record of 
1.024 usee. 

Bicause of the Nyquist criterion, the Doppler sig- 
nal mi st be sampled at least twice per cycle to avoid 
alia: i ng in the reconstructed time ecord. This puts 
an upper limit on the simulated uoppler signal of 0.5 
cycles/sample period. One would prefer to sample the 
signal or the order of 10 times/cyuic, ji°lding an 
effective upper limit of 0.1 cycles/ sample oeriod. 
Furthermore, in order to use discrete Fourier transform 
to perform linear convolutions without aliasing, the saf’ 
of lengths of the nonzero portions of the Doppler burst 
and the nonzero portion of the impulse response of the. 
filter must be less than the length of the tia* rpcoi d 
(6, p. 111). This constraint places a lower limit oi 
the frequencies to be studied using this simulation. 

A ratio of upper to lower cutoff frequencies of ’ 
to 1 was considered general enough to serve the purpos.s 
of this paper. The upper cutoff frequency of the fil- 
ters was chosen as 0.08 cycles/sample period to fall 
well within the upper frequency limit. The lower cutoff 
frequency was therefore 0.04 cycles/sample period. With 
the two cutoff frequencies specified, the impulse 
response of the filter (equal to the inverse discrete 
Fourier transform of the frequency responses of the 
filters) could be calculated. It was found that for the 
four filters, these functions approach zero at or about 
256 sample periods. The nonzero portion of the Doppler 
burst could then extend to 750 sample periods and still 
meet the convolution constraint. This leads to a lower 
frequency limit on the Doppler burst (for Nf = 8) of 
fp = Nf/length = 8/750 = 0.011 cycles/sample period. 
Since the frequencies of interest fall well within the 
range of 0.011 to 0.1 cycles/sample period, the fre- 
quency limits imposed by the Nyquist criterion and the 
discrete Fourier transform are satisfied. 

Counter-Processor Simulation 

The filtered signal, either noise-free or noise- 
added, is then applied to the counter simulation pro- 
gram. The algorithm for the processor is bas’d on a 
commercially available signal processor used at NASA 
Lewis for laser Doppler Anemometry (LDA) measurements 
(7). The processor model performs two primary 
functions - input conditioning and timing. 

Input conditioning can be viewed as defining the 
envelopes used in making the measurement. These enve- 
lopes are the burst envelope, the envelope of the number 


of cycles used in calculation the average (the N f 
envelope), and the comparison envelope. The burst enve- 
lope detector determines the beginning and end of a 
Doppler burst using the thresholds TH1 and TH2 as 
shown in Fig. 4. The lower counter threshold TH1 is 
set to K times the background noise level, where 
is an arbitrary constant. The high threshold TH2 is 
usually set to five times THl. The beginning of the 
burst is defined as the first Doppler cycle which 
crosses TH2. The burst continues until a cycle crosses 
THl but not TH2. In Fig. 4, 'he burst envelope encom- 
passes six valid cycles. The N c envelope starts at 
the first valid cycle and continues until N c cycles 
are found. This is shown in Fig. 4 for Nq = 4. If 
the burst ends before N c cycles are found, the burst 
is ignored by the processor. The comparison envelope is 
defined in the same way as the N c envelope tut with 
fewer cycles. The ratio of the number of cycles in the 
comparison envelope to the number of cycles in the N c 
envelope (called the comparison ratio r) is generally 
set to 1/2 for N c = 2 and N c = 4 ,<nd to 5 '8 for 
N c = 3. Figure 4 shows the comparison envelope 
applicable for N c » 4. 

The timing portion of the model uses the negative- 
going zero-crossings of the valid cycles (circled loca- 
tions in Fig. 4) to determine the Doppler frequency. 
These zero-crcssings are found using linear interpola- 
tion between the time domain samples that bracket the 
zero-crossing. In the actual processor, Schmitt trig- 
gers are used to find the sero-crossings. The Doppler 
frequency is then calculated as 


where T^c is the time interval between the first and 
last zero-crossings within the N ? envelope. 

To determine if the Doppler frequency is changing 
drastically from cycle to cycle within the burst, a com- 
parison test is made. Tr is defined as the time 
interval between the first and last zero-crossing in the 
comparison envelope. In the actual processor, the meas- 
urement passes the comparison test if 


rT 


Nc 


- K < T N < rT Nc 


♦ K 


(5) 


where K is a constant number of clock counts set from 
the front panel. For the simulation program, the com- 
parison test uses the relationship 

l rT Nc " T n| < com P arison accuracy (6) 


where the comparison accuracy used is generally set to 
5 percent. The simulation comparison test is percentage 
based, while the actual processor measures the differ- 
ence in counts. What this implies is that in the real 
processor, as the Doppler frequency increases, the tol- 
erance on the comparison validation increases because 
the limits are a fixed number of clock counts. In the 
digital simulation, the comparison validation is given 
as a percentage of T|y c . The comparison was performed 
in this manner so as to generalize the results of this 
’aper. To relate the percentage to the actual processor 
switch setting, the Doppler frequency must be specified. 

Signal to Noise Ratio (SNR) 

Since the noise output of the photomultiplier is 
assumed to be white, one cannot speak of SNR until the 
signal has been filtered. Herein, we define the peak 
SNR as the ratio of the mean signal power at the filter 
output to the filtered mean noise power where both are 
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averaged over one cycle, with the exponential term 
appearing in Eq. (1) set to its peak value, unity. 
Therefore, the mean signal power can be written as 

P s - 0.5(A,V) 2 |H(f)| 2 (7) 

kith the assumption of shot effect induced white noise, 
the power spectral density (at a given time) is equal 
to the expected photoelectron rate. The noise power is 

P n . 2<x(t)> N bw (8) 

where the noise bandwidth Nfc w is 

N bw * "? /" l H ( f )| 2(lf < 9 > 

G 4 

and G is defined as the maximum value of |H(f)|. 

Note that the rms noise (ensemble average at a given 
time) is 


o 



’he noise power averaged over one cycle is thus 

P„ - 2N. (A„ ♦ A \ 

n bw\ dc sj 

and the peak SNR is 


SNR 


(A s V) 2 |H(f)| 2 

4N bw( A dc 


( 10 ) 


(ID 


( 12 ) 


RESULTS 

The output of the filter was determined by finding 
the inverse Fourier transform of the product of the 
Fourier transform of the input signal and the Fourier 
transform of the filter. The noise-free signals were 
evaluated primarily to determine the ability of the 
filters to attenuate the pedestal without atte'-Ling the 
Doppler frequency. As shown in Fig. 2, pedestal leakage 
is a strong function of Nf. However, it is not clear 
whether the pedestal leakage is a problem in determining 
the Doppler frequency. 

A more obvious source of error in measurements is 
the number of cycles used in calculating the Doppler 
frequency. Figure 5 depicts the percent error as a 
function of frequency of the average of four fringes 
(i.e., N c = 4) for Nf = 4, 8, and 12. It is appar- 
ent that one should not choose Nr to be equal to 
N,., even if the flow is parallel to the fringe normals. 
Either the amplitude is not large enough to allow the 
four cycles to be counted (as is the case with the 
Bessel filter. Fig. 5(b)), or the error in the later 
cycles is large enough to induce relatively urge errors 
in the average (as happens with the other filters). 
Otherwise, the general trend is that the larger the num- 
ber of fringes, the more accurate the measurement will 
be. Other considerations, such as the size of the meas- 
urement volume, force the number of fringes to be as 
small as possible. It was our experience that if Nf 
is chosen to be roughly twice N c , the measurement 
accuracy is sufficient. For purposes of this report, 
it was determined that the combination of Nf ■= 8 
ano N c = 4 was adequate. These values will be used 
throughout the remainder of the report. 


Another observation can be made from Fig. 5. The 
Butterworth filter frequency response exhibits a marked 
slope as a function of frequency. This implies that the 
measured frequency will tend to be biased toward the 
center of the filter passband. It is desirable that the 
slope of the error curve be very nearly flat over a wide 
range in the passband. The Bessel and both Chebyshev 
filters more nearly approach this condition. However, 
the slopes of the error curves in the Chebyshev filters 
become pronounced at the passband edges. 

Table I shows the normalized frequency for indi- 
vidual cycles of the filtered bursts as a function of 
both frequency and filter. The amplitude was increased 
in each case until 10 cycles were accepted by the 
counter processor. In general, the later cycles result 
in a lower frequency than the earlier cycles, and the 
bursts of higher input frequency are biased toward the 
lower frequencies. The exceptions to these general 
trenos are the two Chebyshev filters. At the higher 
frequencies, the later cycles are of higher frequency 
than the earlier cycles. This does not appear in the 
averages because bursts of high enough amplitude to be 
processed first cross TH2 (see Fig. 4) at cycle 4 or 
earlier. Those bursts which cross the upper threshold 
at cycle 5 or higher will not have 4 consecutive cycles 
at greater amplitude than TH2. 

As noise is added to the Doppler signal, a number 
of observations can be made. Figure 6 is the comparison 
of normally distributed random amplitude signals (mean 
value equals 10 counts/sample period, standard deviation 
equals 10 counts/sample period), both noise-free and 
noise-added, as a function of frequency. For the noise- 
free case, the Bessel filter and the 0.1-dB Chebyshev 
filter appear less sensitive to amplitude fluctuations 
than the other two. But, again, for the four filters, 
tne general trend is downward as the frequency 
increases. With Poisson-noise added, the signal accu- 
racy deteriorates, as expected. The “bumps" in the 
Bessel filter frequency response occur because fewer 
Bursts are accepted by the pro essor for this filter 
than for the other three. Sir .e the Bessel filter 
attenuates the input signal, ignals of higher input 
amplitude are necessary to C' ss the upper threshold. 

The data rate using the Bessel niter will be necessar- 
ily lower than for tne other filters. The 0.1-dB 
Chebyshev filter exhibits the flattest profile through- 
out the frequency range. 

Another consideration in the use of laser anemo- 
meters is the measurement of turbulence intensity. 

Figure 7 shows the standard deviation of the measure- 
ments as a function of frequency for the four filters 
for random amplitude signal with noise added. Tne four 
filters are clustered around the 1- to 1. 5-percent 
range. This implies that turbulence intensities of this 
magnitude cannot be measured accurately using the stand- 
ard techniques. Also, measurements of mean velocity 
will be biased near the band edges because of the vari- 
ation of the measurement rate with frequency caused by 
filter attenuation. For these measurements, the Bessel 
filter exhibits the highest standard deviation at the 
higher frequencies. It is desirable that this curve be 
nearly flat to avoid errors in turbulence intensi*y as 
a function of frequency. The flattest profile can be 
seen to be the Bjtterworth filter, especially in the 
center of the passband, although the differences between 
these filters in this respect are not especially 
significant. 

Figure 8 shows the SNR of the four filters as a 
function of frequency for fixed-amplitude, noise added 
input signals. The Bessel filter exhibits the least 
variation in SNR across the passband, but the SNR is 
significantly lower than for the ether filters. The 
0.1-dB Chebyshev filter also has a fairly flat profile, 


4 



but at higher magnitude than the Bessel. Both the 
Butterworth and the 1-dB Chebyshev filters exMbit a 
large variation of SNR as a function of frequency. 

This implies that the input signals will be biased more 
strongly by these filters than by the Bessel and the 

0.1-dB Chebyshev filters. This observation is borne 
out by the graphs of Fig. 6. 

CONCLUSIONS 

From this admittedly limited study, several impor- 
tant conclusions can oe reached. First, the error in 
the signal can be viewed to be partially a function of 
the ratio bet^aen Nf and N c . As a rule of thumb, 
if Nf/N c « the errors caused by the large vari- 
ation in cycle frequencies at the end of the burst will 
be minimized. Second, for the four filters studied, 
there is generally strong signal biasing near the band 
edges. This can be avoided by defining "effective" 
cutoff frequencies inside the actual cutoff frequencies 
of the filters. A 10-percent s *• i f t upwards for the low 
edge, and downwards for the high edge, should be suffi- 
cient to limit these errors. Third, biasing errors 
caused by the filters can be linked to the variations 
in SNK across the passband. A filter with a flat SNR 
profile is then more desirable than a filter with a 
large variation in SNR as a function of frequency. 
Fourth, variations in standard deviations, which can 
cause errors in turbulence intensity measurements, can 
be interpreted as a function of tne magnitude of the SNR 
of the filtered signal. F or this reas:n, and for the 
obvious reason that a higi.er SNR yields a higher da. a 
rate, a filter should be chosen which allows as high an 
SNR as possible. With these four points in mind, the 
Chebyshev filter with the 0.1-dB ripple is the best 
filter of the four studied. 

It should be noted that this study is by no means 
exhaustive. No effort was made to determine the errors 


Induced by these filters on the measurement of turbulent 
flow or of flows for which the mean velocity vector 
crosses the measurement volume at an angle to the fringe 
normals. Furthermore, there exist a large number of 
filters not evaluated here as well as higher order fil- 
ters of the types examined. A more extensive study of 
these systematic errors is required to minimize the 
inaccuracies of the measurements made by fringe-t.vpe 
laser anemometry systems. 
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TABLE 1. - NORMALIZED FREQUENCY FOR INDIVIDUAL CYCLES OF FILTERED SIGNAL 


BE * Bessel 
BU * Buttermorth 
CH * Chebyshev 

Number fringes (Nf) * 8 


BANDPASS FILTER CUTOFF FREQUENCIES 
FI = >040 cycles/sample period 
Fu = .080 cycles/sample period 


Fd filter CYCLE 


cycles/ 


sample 

period 



1 


2 


3 


4 


5 


6 


7 


8 


9 


’0 

.040 

BE 


1. 

008 

1, 

007 

1. 

006 

1. 

005 

1 • 

004 

1. 

003 

1. 

001 

l , 

000 

0. 

998 

0. 

997 


BU 


1. 

017 

1. 

,017 

1. 

016 

1. 

016 

1. 

015 

1. 

014 

1. 

012 

1. 

010 

1. 

008 

1. 

004 


CH- 

ldB 

1. 

026 

1. 

,026 

1. 

026 

1 . 

025 

1. 

024 

1. 

021 

1. 

,016 

1. 

014 

1. 

009 

1. 

,00* 


CH- 

. ldB 

1. 

017 

1. 

,016 

1. 

014 

1. 

012 

1. 

008 

1. 

003 

0. 

996 

0. 

986 

0. 

970 

0. 

948 

.045 

BE 


1. 

006 

1. 

,005 

1. 

004 

1. 

003 

1. 

002 

1. 

001 

1. 

000 

0. 

999 

0. 

997 

0. 

,995 


BU 


1. 

012 

1. 

,011 

1 • 

010 

1. 

009 

1. 

007 

1. 

005 

1. 

,002 

0. 

997 

0. 

992 

0. 

,984 


CH- 

ldB 

1. 

016 

1. 

,015 

1. 

013 

1. 

009 

1, 

C04 

0. 

,996 

0. 

,963 

0. 

,9*3 

0. 

930 

0. 

,891 


CH- 

.ldB 

1. 

Oil 

1, 

,009 

1. 

007 

1. 

,004 

0. 

,999 

0 

994 

0. 

,987 

0. 

,977 

0. 

963 

0. 

,941 

.050 

BE 


1. 

,004 

1, 

,004 

1. 

,003 

1, 

,002 

1 . 

,001 

1. 

000 

0, 

,999 

0. 

,998 

0. 

,996 
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Fig. 1 Digital time records and frequency spectrum simulations of Doppler burst (fQ * 0.060; A(j c ■ 10 
counts /sample period; A $ * 10 counts /sample period; V -1.0; N f = 8; t Q “ 512 ; v>* 0.23475 rad.) 

la) Noise-free time record (b) Noise-free frequency spectrum 

(c) Poisson noise-added time record (d) Noise-added frequency spectrum 

(e) Time record of noise-added signal (f) Frequency spectrum of noise-added 

filtered by Chebyshev filter (0.1 dB) filtered signal 

(F| - 0. 04 cycles /sample period ; 

F u • 0.08 cycles /sample period) 
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Fig. 2 Frequency spectrum ol Doppler burst 
shewing effect of Nf on pedestal leakage 
(fp ■ 0.08 cycles /sample period) 

(a) Unfiltered frequency spectrum 
with Nf ■ 4 

(b) Unfiltered frequency spectrum 
with N f ■ 12 

(cl Signal from (a) filtered with 
Bessel filter (F| • 0.04 cycles/ 
sample period; F u -0. OB cycles/ 
sample period! 

(di Signal from (b) Mitered with 
Bessel filter (Fj * 0.04 cycles/ 
sample period; Fy ■ 0. OB cycles/ 
sample period) 






Fig. 4 Measurement of Doppler frequency by 
counter-processor (processor measures elapsed 
time for zero crossinqs ;f Nj- cycles that cross 
upper threshold (TH2I. Shown are burst envelope, 
N c envelope, and comparison envelope. ) 




PERCCNT ERROR 





Fig. 6 Percent error as function of frequency 
of the averages of four fringes (Nf • 8, N c *4) 
for normally distributed random amplitude 
signals 

(;' Butterworth filter 
<b) Bessel filter 

ic) Chebyshev filter (1-dB ripple 
(d) Chebyshev filter (0. 1-dB ripple) 
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Fig. 7 Standard deviation of the measurement of frequency as function 
of frequency of the average of .'our fringes (N f ■ 8, N c * 4) for nor- 
mally distributed random amplitude signals with noise added 



Fig. 8 Signal to noise ratio as function of frequency for A s - 15 
counts /sample period, A^ - 10 counts /sample period, noise- 
added signals 


1. Report No. 2. Government Accession No. 

NASA TM-87047 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Filter Induced Errors in Laser Anemometer 
Measurements Using Counter Processors 

5. Report Date 

6. Performing Organization Code 

533-04-1A 

7. Authors) 

Lawrence G. Oberle and Richard G. Seasholtz 

8. Performing Organization Report No. 

E-2609 

10. Work Unit No. 

9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


Prepared for the International Laser Anemometry Symposium, ASME Winter Annual 
Meeting, Miami, Florida, November 17-21, 1985. 


16 Abstract 

Previous simulations of laser Doppler anemometer (LDA) systems have focused pri- 
marily on noise studies or biasing errors. Another possible source of error is 
due to the choice of filter types and filter cutoff frequencies. In general, 
before it is applied to the counter portion of the signal processor, a Doppler 
burst is filtered to remove the pedestal and to reduce noise in the frequency 
bands outside the region in which the signal occurs. Filtering, however, can 
introduce errors into the measurement of the frequency of the input signal which 
leads to inaccurate results. With the results of this paper it is possible to 
evaluate errors caused by signal filtering in an LDA counter-processor data 
acquisition system and to choose filters for a specific application which will 
reduce these errors. 


17 Key Words (Suggested by Authors)) 

Laser anemometry 

18. Distribution Statement 

Unclassified - unlimited 
STAR Category 35 

19. Security Classil. (ol this report) 

Unclassified 

20. Security Classlf. (of this page) 

Unclassified 

21. No. of pages 

22. Price* 


‘For sale by the National Technical Information Service, Springfield, Virginia 22161 



















